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Abstract 

We assume a community whose members adopt one of two opinions. Each member 
appears as an inflexible, or as a non-contrarian or contrarian floater. An inflexible 
adheres to its opinion, whereas a floater may change into a floater of the alternative 
opinion. The occurrence of this change is governed by the local majority rule: members 
meet in groups of a fixed size, and a floater changes provided its opinion has the 
minority in its group. A non-contrarian floater keeps the opinion it adopts under 
the local majority rule, whereas a contrarian floater adopts the alternative opinion. 
We determine the dynamics generated by repeated opinion changes, and study its 
dependence on the densities of inflexibles and the frequencies of contrarians among the 
floaters. We restrict ourselves to groups of size at most 3. 

Keywords: Sociomat hematics, sociophysics, opinions dynamics, local majority rule, con- 
trarian behaviour, floating behaviour 

NOTE concerning the current version of the text: due to the large sizes of some 
figures, the accompanying captions to these figures do not fit together with the 
figures on the same page. Therefore, these captions are placed in the main text, 
instead of in the sequence of Figures at the end of the text. 

1 Introduction 

Within the growing field of sociophysics (see [T] for the defining paper, and [3]-j6] for an 
impression of the recent state of the art), two- and three states Galam opinion models and 
their modifications (j3j, [6]) play a guiding role in analyzing the process of opinion spreading 
in communities. These models are centered around the local majority rule (l.m.r.), which 
is applied either in a deterministic or a probabilistic way. In the basic deterministic case, 
supporters of the (two or three) opinions present in a community are randomly distributed 
over all possible groups of a fixed size. Within each group the members adopt the opinion 
that has the majority in that group, after which all group members are recollected again. 
In case there is no majority in a group, its members stick to their own opinion (i.e., neutral 



3 



treatment; the probabilistic treatment in case of a tie assigns opinions to the group mem- 
bers according to a certain probabihty distribution). Repeated apphcation of this principle 
generates what is caUed randomly localized dynamics with a local majority rule ([6j). In the 
basic probabilistic case, the community members are divided among groups of various sizes 
according to some probabihty distribution, and within each group aU members adopt one of 
the possible opinions with a certain probability. 

In the basic deterministic two states opinion model, fast dynamics occurs in which the 
opinion that originally has the majority eventually will obtain complete presence at the cost 
of the alternative opinion. In the probabilistic two states opinion model, the final outcome 
depends on the probability distributions for group sizes and local adaptation. Eventually 
the state of the community can be either a polarization on the initial majority or minority, 
or a perfect consensus on both opinions (see [3j , which unifies basic probabilistic two states 
opinion models). In [6] a three states opinion model is introduced in which the community 
members are randomly distributed over groups of size 3. Within each group the l.m.r. is 
applied, with the additional rule that in case of a tie all members of the group adopt one of 
the three opinions according to some probability distribution. It is shown that the dynamics 
leads to fast polarization on one of the three opinions, which may be an opinion that initially 
has a minor presence in the community. 

To add more realism to opinion models, in [2j the basic deterministic two states Galam 
opinion model is extended by the introduction of so-called contrarians. A contrarian is a 
community member who, instead of the opinion it should adopt under the l.m.r., switches 
to the alternative opinion. Depending on the density of contrarians as well as on group size, 
their presence either leads to a stabilization of the opinion dynamics in which one opinion (the 
one with the lower density of contrarians) dominates the other, to an equilibrium in which 
neither opinion dominates (in case both opinions have equal densities of contrarians), or (in 
case of relatively large densities of contrarians for both opinions) to a dynamics in which 
the dominating opinion constantly alternates between the two opinions. The incorporation 
of contrarians in opinion dynamics models is a step towards an explanation of the "hung 
elections" outcome in several recent voting events (the U.S.A. presidential elections in 2000, 
the German elections in 2002 and 2005, as well as the 2006 Italian elections, see [2j, [5j), in 
which one strategy achieved only a very narrow gain at the cost of another strategy. The 
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origin of contrarian behaviour as well as its implications have been the focus of various other 
recent studies, see [8], [H], [I2], [13]. 

In a recent paper ([4]) the basic deterministic two states Galam model has been modified 
by introducing opinion supporters that express what in politics (and other games) is called 
inflexible behavior. An inflexible community member is a supporter that under all conditions 
sticks to its opinion. Under this terminology supporters that switch opinion when in the 
minority then classify as floaters^ and we shall use this distinction in what follows. In [4] 
the effect of inflexible behavior on opinion dynamics is studied for the case that opinion 
supporters repeatedly meet in groups of flxed size 3. It is shown that a small density of 
inflexibles for only one of the two opinions allows for the existence of two local attractors. 
One of these local attractors is a mixed one, on which both opinions are present and on which 
the opinion that is supported by inflexibles is a minority. The other attractor is a single 
state attractor, on which the opinion that is supported by inflexibles has complete majority, 
i.e., its density equals 1, the other opinion being absent. Due to the presence of these two 
attractors, the outcome of the opinion dynamics thus depends on the initial condition, the 
basin of attraction for the mixed local attractor being relatively small compared to that for 
the single state attractor. If the density of inflexibles is sufliciently large (approximately 
17%), the mixed attractor disappears and the single state attractor becomes global. In 
case both opinions have small and equal densities of inflexibles there are two mixed local 
attractors. These two attractors are symmetrically situated with regard to a separator on 
which both opinions are present with density 0.5. A change in the density of inflexibles for 
one of the opinions breaks this symmetry, and a sufliciently large increase may lead to a 
global attractor on which the opinion with the larger density of inflexibles has the majority, 
see [4j. Recently, the effect of inflexibles and floaters on opinion dynamics has also been 
studied in [9j and fTO\ . 

In this paper we combine the approaches presented in [2j and [4J, by allowing for groups 
composed of inflexibles as well as contrarian and non-contrarian opinion supporters. For 
clarity we restrict ourselves to groups of flxed size at most 3. For both opinions we assume 
flxed densities for the inflexibles. Also, we consider the contrarians to be part of the floaters, 
i.e., in a given group the contrarians flrst determine their opinion according to the l.m.r., and 
subsequently change to become a floater (not necessarily a contrarian) for the alternative 
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opinion (which thus may be the opinion that the contrarian initially was supporting). The 
presence of contrarians for each opinion is quantitatively expressed as a fixed fraction of the 
density of fioaters of the respective opinion. In case of a tie in groups of size 2 we apply the 
neutral treatment. After an opinion update, all supporters for both opinions are recollected 
and then are redistributed again, either as an infiexible or as a non-contrarian or contrarian 
floater, according to the fixed densities for infiexibles and the fixed fractions of contrarians 
for the two opinions. We study qualitative characteristics of the opinion dynamics gener- 
ated by repeated updates. In particular we study changes in the number of equilibria, and 
changes from monotone to alternating dynamics, due to changes in parameter combinations. 
A detailed mathematical treatment of the generated opinion dynamics, for groups of arbi- 
trary size, shall be given in a forthcoming paper (^). 

Notation 

We denote the two opinions by A and B. The densities of infiexibles for the A and B opinion 
are denoted by a a and as respectively, with < ^ 1 as well as < < 1, and in addi- 
tion < aA + oiB ^ 1- Since the roles of the A and B opinion are interchangeable in deriving 
the opinion dynamics, we may without loss of generality assume that < ^ 0.5, and we 
shall do so in what follows. The fraction of contrarians among the A fioaters is denoted by 
7^, and 7^ denotes the fraction of contrarians among the B fioaters, with both < 7^ < 1 
and < 7b < 1. The size of the groups in which opinion supporters meet is denoted by L. 
The density of the A opinion at time t = 0, 1, 2, • • • (or after t updates) shall be denoted as 
Pt. Note that for given a a and as the density pt necessarily lies in the interval [a^, 1 — as] 
(independent of L, 7^ or 7^). With fL;aA,aBnA,iB denote the function that determines 
the density of the A opinion after application of the l.m.r. followed by the switch of the 
contrarians. Thus, pt+i = fL^aA^aBHAnBiPt)- Setting 7a = 7b = 0, Pt+i = fL;aA,aB;0,o{Pt) 
then gives the density obtained from pt when the l.m.r. is applied without being followed 
by the switch of the contrarians. In the Appendix tables are given, presenting all possible 
group compositions in terms of infiexibles and non-contrarian and contrarian fioaters for 
group sizes L = 1 to 3, together with the effects of the l.m.r. and the opinion changes of 
contrarians. It is assumed that the community is sufficiently large to allow for the derivation 
of the density of each possible group composition in the ensemble of all groups of a fixed size 
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from the densities in the community of the constituents of a group. From these tables the 

expressions for /L;aA,«s;7A,7s obtained. 
> 

With fL-oLA,oLB]iA,iB denote the dynamics generated by repeated apphcation of fL^aA^aBHAHB 

in subsequent timesteps. Furthermore, Ph-^aA^o-B^iA^iB denotes an asymptotically stable equi- 

> 

librium for fh-cxA^oLB^iA^iB^ Pl-ola,olb]ia,ib ^^^^^^ ^o an asymptotically stable periodic point. 
We now turn to the treatment of the opinion dynamics for group sizes L = 1, 2 and 3. 



2 L = 1 

The case L = 1 resembles a community in which each member is unaffected by other com- 
munity members in determining its opinion, and the only changes in opinion come from the 
contrarians. The contributions to the A density after application of the local majority rule 
is obtained from the second column in Table 1 in Appendix |6.1[ This column obviously is 
equal to the first one, since in groups of size 1 local majority is automatically obtained, but 
is without eflFect on the opinion densities. These contributions are: aA for the A inflexibles, 
and p — oiA for the (non-contrarian and contrarian) A fioaters. Their sum is and we obtain 
that 

Pt+i = /i;aA,c.s;0,obt) = Pt] (1) 

consequently, eachp G [a a-, I — olb] is a neutrally stable equilibrium for the opinion dynamics. 
In case only (non-contrarian and contrarian) fioaters are involved both aA and as are equal 
to 0, and we restrict ourselves to the contributions from the second, third, fifth and sixth line 
in the table. Since the l.m.r. leaves each group of size 1 unaflFected, a switch by a contrarian 
in this case necessarily implies a change to the opinion it initially does not support. Thus, 
here also a contribution to the A density comes from the group that initially consists of only 
B contrarians, as these will turn into A fioaters. In this case we obtain for the contribution 
to the A density: 

Pt^l = /l;0,0;7A,7sbt) = (1 " 7A)Pt + 7^(1 " Pt) = 7B + (l " (TA + 7^))^^. (2) 

The effect of both infiexibles and non-contrarian as well as contrarian fioaters is obtained 
by adding all the expressions in the last column: the contributions aA due to the invariant 
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density of A inflexibles, (1 — 7a) (Pt — from the non-contrarian A floaters, and 7^(1 — 
— Pt) from the contrarian B floaters. This yields: 



Pt+i = hv 



\OLA,OLB]lA,lB 



iPt) 



ttA + (1 - 7A)(Pt - c^a) + 7b(1 - ttfi - Pt) 



(3) 



«a7a + (1 - 0:3)73 + 1 - (7a + 7b) )Pt- 



It follows that if 7a + 7b > 0, then 



P = 



oaIa + (1 - «b)7b 
lA + 7b 



(4) 



is the unique equihbrium for the opinion dynamics fi-aA,(^B-,iA,iB- ^o its hnearity as a 

function of p^, expression ^ imphes that the dynamical characteristics of this equilibrium 
are governed solely by the frequencies of the contrarians. The equilibrium is asymptotically 
stable if and only if < 7^ + 7b < 2. For < 7^ + 7b < 1 the equilibrium is approached 
monotonically, with an increase in the A density if and only if its initial value is less then 
the equilibrium value. For 7a + 7b = 1, the function fi;aA,aB;7A,7B constant and equals 
<^a7a + (1 ~ oib)ib] the opinion dynamics then reaches its equilibrium in one iteration. For 
1 < 7a + 7b < 2, the equilibrium is approached alternately. For 7a + 7b = 2, i.e., both 
7a = 1 and 7b = 1, the equilibrium equals 0.5(1 + q^a — <^b) and is neutrally stable; each 
p G [q^a, 1 — <^b] different from 0.5(l + aA — <^b) generates a neutrally stable cycle of length 
2. 

On the equilibrium, the A opinion has the majority if and only if the inequality 



holds. Thus, for an opinion to achieve the majority it is required that it is being supported by 
a sufliciently large density of inflexibles, and/or a sufliciently small frequency of contrarians 
among the floaters. 

Given densities a a and q^b of inflexibles for the two opinions, a change in the frequencies of 

contrarians from into small values 7a and 7b causes the bifurcation from a collection of 

> > 

neutrally stable equilibria for /i;a^,a^;o,o into a unique stable equilibrium for fi-oLA^oLB^iA^iB- 



(0.5 - aA)lA < (0.5 - aB)7B 



(5) 
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The opinion which has the majority on this equihbrium is determined by inequahty ([5]). 
In case = = a, the opinion with the smaller frequency of contrarians obtains the 

majority. Conversely, given different frequencies 7^ and 7^ of contrarian floaters for the 

> 

two opinions, in the absence of inflexibles the dynamics /i;o,o;7a,7s has p = ^^f^^ as its 
unique stable equilibrium, on which the opinion with the smaller frequency of contrarians 
has the majority. Fixing sufliciently small densities a a and of both opinions as inflexibles, 
this equilibrium slightly shifts but leaves the majority unaltered. In case 7a = 7b, in the 
absence of inflexibles the equilibrium p equals 0.5, and the introduction of small densities of 
inflexibles for both opinions changes this equilibrium into one on which the opinion with the 
larger density of inflexibles takes the majority. Figure [T] illustrates these conclusions. 

[Figure 1 about here.] 

Figure [2] gives a qualitative overview of the outcomes of the possible opinion dynamics 

> 

[Figure 2 about here.] 

Figure [2| An overview of the opinion dynamics of fi;aA,aBnA,iB^ values aA and as as 
indicated, and with ja and jb in each pane on the horizontal and vertical axis, respectively, 
both ranging between and 1. In each pane the line with negative slope 7a + 7b = 1 
is drawn, and possibly an additional dashed line of positive or zero slope. On the line 
7a + 7b = 1 the function fi-aA.dB^iA^iB i^ constant, and the corresponding values of 7^ 
and 7^ separate between monotone and alternating dynamics, with the monotone dynamics 
occurring if < 7^ + 7b < 1, i-e., below the line. The dashed line, if present, gives the 
values (71,72) 7^ (0,0) for which the equilibrium of the opinion dynamics equals 0.5, and 
is determined by the expression [aA — 0.5)7a — {olb — 0.5)7^ = 0. Opinion A obtains the 
majority if (and only if) [aA — 0.5)7a — {olb — 0.5)7^ > holds, i.e., if as < 0.5 and (7^, 7b) 
lies above the dashed line. The panes for values (a^, <^b) for which aA + olb = ^ represent 
degenerate cases, in the sense that these values allow for only one density p = aA for the A 
opinion to occur. 
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3 L = 2 



In groups of size 2 the number of members that support the A oi B opinion may be equal, 
in which case a tie occurs. We shaU deal with the neural treatment in case of a tie, in which 
each supporter keeps its own opinion. 
Appendix |6.2| shows the table for this case. We obtain 

Pt^l = /2;aA,«s;0,obt) = Pt, (6) 

which is obvious, since in groups of size 2 no majorities can occur, and, in case of a tie, the 
neutral application of the local majority rule does not have any effect. Incorporating the 
effect of non-contrarian as well as contrarina floaters, the table yields that 

Pt+l = f2;aA,aB;7A,7B{Pt) = ^^AIA + (l " (^b)1B + " (TA + 7^))^^. (7) 

Thus, for groups of size 2 the effect of the neutral application of the local majority rule and 
the contrarians is the same as for groups of size 1. 

4 L = 3 

Group size 3 is the smallest value of L for which the local majority rule becomes effective 
due to possible group compositions in which a majority of one of the two opinions occurs. 
As a consequence, the generated dynamics allows for features different from those for group 
sizes 1 and 2. Careful bookkeeping based on the table in Appendix |6.3| yields that 

iPt) 

= q;a(1 - 7a) + (1 - «b)7b - (2«^(1 - 2ja) - + iB^Pt 

+ (3 + aA{l - 27a) - - 275) - 47^ - 275)^? _ 2(1 - 7^ - 75]^? 

= + q;a(1 - 7a) + (1 - «i3)7B ~ (1 + 2«a(1 - 27a) - 7a + iB^Pt 

+ (3 + aA{l - 27a) - oib{1 - 27b) - 47a - 27b)p? - 2^1 - 7^ - 7b)p?- 
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For clarity we start the analysis of the generated opinion dynamics with the symmetric case 
of equal densities of inflexibles and equal fractions of contrarians for both opinions. 

4.1 The fully symmetric case: aA = c^b sind 7^ = 75 

Taking aA = c^b = and 7a = 7b = 7, we obtain that 

Pt+l = fs;a,a;^M = + (l " 2^,) (7 + a{l - 2j) - (1 - 2j)pt + (1 - 2j)p^^ (9) 

As an illustration to expression |9| Figure [3] shows a collection of graphs of fs-a,a;-f,-f as function 
of for a = 0.1 and several values of 7. 

[Figure 3 about here.] 

From expression ^ the analysis of the generated opinion dynamics is straightforward. We 
give an overview. 

> 

Symmetry considerations imply that the dynamics fs-a,a;-f,-f has p = 0.5 as an equilibrium, 
for any choice of a G [0, 0.5] and 7 G [0, 1]. In addition to parameter combinations a and 7 
for which this equilibrium is unique and stable, there are combinations which allow for an 
unstable repelling equilibrium p = 0.5 in combination with two other, asymptotically stable, 
equilibria, or with two asymptotically stable periodic points of minimal period 2. Details for 
these possibilities to appear are derived in Appendix |6.4[ here we confine ourselves to the 
outcome. 

Let the critical curves C3 and Cs be defined as follows: 

C3 = {(«,7) G [0,0.5] X [0, 1] : (3 - 4«)(1 - 2j) = 2}, (10) 

and 

Cs = {{a, 7) G [0, 0.5] X [0, 1] : (3 - 4a)(l - 27) = -2}. (11) 

Figure [i] shows the curves C3 and Cs in the (a, 7)-parameter space. On C3 the derivative 
/3;a,a;7,7(0-5) ^quals 1, whcrcas on Cs this derivative equals -1. The two corner areas in 
Figure [4] enclosed by either Cs or Cs are the regions of parameter combinations for which 0.5 is 
unstable; outside these regions (including the curves) 0.5 is the unique asymptotically stable 
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equilibrium for fs-a.a-^-r.-ri independent of the initial condition. The lower left corner region is 

y 

the area for which the dynamics /3;a,a;7,7 has two asymptotically stable equilibria Ps-a.a^j.j- 
Given parameter combinations (a, 7) in this region, the opinion dynamics eventually will 
stabilize on an equilibrium on which the opinion with the initial majority will have maintained 
its majority. In case (0^,7) 7^ (0,0), this equilibrium is mixed; if neither inflexibles nor 
contrarians are present for both opinions, i.e. (a, 7) = (0, 0), the equilibrium is a single state 
attractor with only one opinion present. These results generalize those obtained in [4J for 
the case of equal densities of inflexibles and no contrarians for both opinions. For parameter 
combinations in the upper left corner region in Figure [4| the dynamics has two attracting 
periodic points of period 2. Here an initial majority does not guarantee the eventual majority, 
since the dynamics is such that both opinions alternately switch between minority and 
majority. 

[Figure 4 about here.] 

Thus, if both opinions are being supported by equal densities a of inflexibles and equal 

fractions 7 of contrarians among the floaters, for an opinion to obtain the majority it is 

necessary that a as well as 7 are sufliciently small, and that it has the initial majority. 

Also, with increasing a (7), the maximum value of 7 (a) for which a majority is attainable 

decreases. If no inflexibles are present, the fraction of contrarians among the floaters must 

be less than approximately 17% ( %) for a majority to be realizable, and if the fraction 

6 

of contrarians among the floaters equals 0, the density of inflexibles must be less than 25%. 
If in the parameter space a combination (a, 7) approaches from within a corner towards 
one of the two critical curves, then the two additional equilibria or periodic points approach 
towards p = 0.5] a withdrawal in the parameter space results in the opposite movement 
of the additional equilibria or periodic points. It follows that when passing through C3, the 



dynamics f3;a,a;j,j undergoes a supercritical pitchfork bifurcation, and when passing through 
Cs the dynamics undergoes a period doubling bifurcation (flip bifurcation). 

4.2 The general case 

We now return to the general expression g and give an overview of the possible outcomes 
> 

of the dynamics /3;aA,aB;7A,7B- The analytical background is given in Appendix 
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6.5 



We 



distinguish several cases. 

For 7^ and 7^ such that 7a + 7b = 1, the function fs^aA^aBHA^iB quadratic in p. 

> 

The corresponding opinion dynamics fs^aA^^BnA^-iA ^ unique stable equilibrium 
in the interval [a^, 1 - c^b]- For (7a,7b) = (0.5,0.5), the function /3;aA,«s;0.5,o.5 be- 
comes constant and equals f3;aA,aB;0.5,o.5{p) = 0.5(1 + aA — c^b)] it allows for a unique 
stable equilibrium p = 0.5(1 + — a^), on which opinion A has the majority if and 
only if aA > c^B- The following two figures distinguish between parameter combina- 
tions aA^ c^B and 7^ for which the A opinion obtains either the majority or minority 
in equilibrium (Figure [5]), and for which the equilibrium is approached monotonically 
or alternately (Figure [g]). It follows that with increasing value of 7^, the region of 
parameter combinations (a^, as) for which opinion A obtains the majority, decreases. 
In addition, if 7^ < 0.5, the A opinion can obtain the majority for any value of a^, 
provided that as is sufficiently small; if 7^ > 0.5, aA must be sufficiently large and 
as sufficiently small for an A majority to occur. 



[Figure 5 about here.] 
[Figure 6 about here.] 

2. 7a + 7b 7^ 1. 

The expression fs;aA,aBnA,iBip) — P = for determining the equilibria is 

fs;aA,(^B;7A,7B (p) ~ P 



= ^^(1 - 7a) + (1 - c^BhB - (1 + 2q^a(1 - 27a) - 7a + 1b)p+ 

(3 + aA{l - 27a) - anil - 27^) - 47a - 27^)^^ _ 2(^1 _ _ ^^^^3 

= 0. 

The number of solutions is determined by its discriminant, which is denoted by 



(12) 



Z)(q;a, Q^b; 7a, 7b)- The expression for the discriminant is derived in Appendix 6.5 
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here we discuss its implications. 

For parameter combinations (o^a, 7a, 7b) such that i^(o^A, 7a, 7b) > 0, the 
equation f?>-aA,aB;iA,iB{p) — P = ^ has a unique real solution. If D{aA^ oib] 7a, 1b) < 0, 
there are three real solutions. However, these solutions do not necessarily have to be- 
long to the interval [q^a, 1 — <^b] (but if a solution lies in this interval, it clearly is an 
equilibrium for the dynamics fs^aA^aBHAnB)- D{aA^(^B]lAnB) = there are three 
real solutions, of which at least two coincide; if this happens in the interval [aA, 1 — <^b], 
the parameter combination is at a bifurcation point, discriminating between dynamics 
with either a unique equilibrium or three equilibria. If at the bifurcation point exactly 
two of the three solutions coincide, the coinciding solutions form a semistable equilib- 
rium. 

Figure [7| shows a collection of signplots for the discriminant, for values a a and as 
indicated, and with 7a and jb for each signplot between and 1. In addition the 
outcome of the analysis for parameter combinations (o^a,<^b;7a, 1 — 7a) is included. 



as well as the results of the analysis for combinations (a, a; 7, 7). Figure 14 in Ap- 
pendix |6.5| shows an extended version of Figure [7| that includes more values of a a and 
aB- It shows that for both a a and aB larger than 0.3, there is not much qualitative 
change in the signplots of the determinant D{aA^ <^b; 7a, 7b)- 



[Figure 7 about here.] 



Figure [7| A collection of panes, for values aA and aB as indicated, and 7a and 7^ 
for each pane ranging between and 1, with 7a on the horizontal axis and 7^ on the 
vertical axis. In each pane the signplot of the discriminant D(q^a, <^b; 7a, 7b) is shown 
for points (7a, 7b) for which 7a + 7b 7^ 1- White areas represent the parameter com- 
binations with a positive discriminant (i.e., combinations for which the corresponding 
opinion dynamics has a unique equilibrium), and in dark regions the discriminant is 
negative (the corresponding opinion dynamics then has 3 different equilibria, but not 
necessarily in the interval [q^a, 1 — <^b])- In each pane, on the line 7a + 7b = 1 a density 
plot of the equilibrium density P3;aA,aB;7,i-7 is given. An increasing gray level indicates 
an increasing equilibrium value, with graylevel (black) representing a density 0, and 
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gray level 1 (white) representing density 1. Furthermore, in panes for which aA = c^b 

holds, on the line 7^ = 7^ in the black regions (i.e., a negative discriminant) in white 

> 

the points are indicated for which the equilibrium p = 0.5 for fs-^aA^aBHAnB unstable; 
other points on the lines 7^ = 7^ (for aA = c^b) indicate parameter combinations for 
which p = 0.5 is stable (as follows from Figure [i]). 



The discriminant becomes singular for parameter combinations (o^a, 7a, 7b) with 
7a + 7b = 1- In approaching such parameter combinations for which (7a,7b) 7^ 
(0.5,0.5), the value of /^(q^a, <^b; 7a, 7b) goes to — oc. For (7a, 7b) = (0.5,0.5), the 
limit generically equals +oc when this point is approached from the region 7a + 7b < 1; 
the limit equals — oc in case it is approached from the other side, i.e., from the region 
7a + 7b > 1- (In case (0.5, 0.5) is approached along the zero set of D{aA^ olb] 7a, 7b), 
i.e., in each pane in Figures [7| and [14] along the boundary that distinguishes between 
the white and black regions and touches with the line 7a + 7b = 1, the limit clearly 
equals 0.) 

Our further discussion of the opinion dynamics fs^aA^aB-nAHB based on Figures [t] 



and 14, Instead of a detailed analytical treatment, we continue with a number of 



characteristic outcomes of the opinion dynamics. 

A first characteristic that draws attention in Figure [7| is the existence of a wedge-shaped 
region of parameter combinations (q^a, <^b; 7a, 7b) with negative discriminant for sufficiently 
small values of all four parameters. For the cases with both aA = c^b and 7a = 7b within this 
region, we already found the existence of two attracting equilibria, symmetrically positioned 
with respect to a third, unstable equilibrium 0.5. We therefore expect also to find a similar 
pattern of three equilibria in [aA, 1 — <^b] for deviations from such symmetric cases within 
the wedge-shaped region. In [4J it has been derived that this is indeed the case in the absence 
of contrarians, i.e., for parameter combinations for which 7a = 7b = 0, and for a a and as 
suflBciently small. Figure [8| which shows a number of graphs of functions fs^aA^aBHA^iB 
relatively small values q^a, <^b, 7a and 7^, implies the same pattern: in case the determinant 
/^(q^a, <^b; 7a, 7b) is negative, the opinion dynamics has two attracting equilibria that are 
separated by an unstable one. The two attracting equilibria differ with respect to the opinion 
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by which they are dominated. By leaving the wedge-shaped area, a bifurcation in the opinion 
dynamics occurs on its boundary /^(q^a, 7a, 7b) = 0. GenericaUy, when moving from 
inside the wedge-shaped area towards this boundary, the unstable equilibrium and one of 
the two stable equilibria move towards each other, and at the bifurcation point merge (thus 
causing a supercritical saddle-node bifurcation). Once the boundary has been crossed, the 
region of parameters with a positive discriminant is entered, and the dynamics is left with 
one attracting equilibrium. On this equilibrium opinion A dominates if the upper part of 
the boundary is crossed, i.e., when 7^ > 7^; opinion B has the majority when the right- 
hand side of the boundary is passed, on which 7^ > 7b holds. This is also illustrated 
in Figure |8) The occurrence of such a bifurcation may lead to a drastic change in the 
outcome of the opinion dynamics: whereas inside the wedge-shaped region the outcome of 
the opinion dynamics depends on the initial condition, outside the wedge-shaped area the 
opinion dynamics will end on the unique equilibrium, independent of the initial condition. 
At the bifurcation point at the endpoint of the sharp region of the wedge-shaped area a 
supercritical pitchfork bifurcation occurs, in which the three equilibria merge together into 
one attracting equilibrium. 

[Figure 8 about here.] 

Figure [8j Four panes of graphs of functions /3;aA,«B;7A,7s relatively small values a^, 
7a and 7^, with (q^a,<^b) as indicated below each pane. In each of the four panes, 7^ and 
7^ satisfy 7a + 7b = 0.2. The color coding of the graphs is as follows: red corresponds 
to 7a = 0, for green 7a = 0.05, for blue 7a = 0.1, for purple 7a = 0.15, and for magenta 
7a = 0.2. I.e., in the corresponding panes in Figure [7] we traverse the line 7a + 7b = 0.2 
from its upper left point on the 7a = axis to its lower right point on the 7b = axis, thus 
passing through regions with positive, zero as well as negative discriminant. Above each 

pane the densities for opinion A are presented, as obtained by the corresponding opinion 

> 

dynamics f?>-aA,aB]iA,iB^ with initial density p = 0.5. 

The white regions in Figure [7] (and also in Figure [14]) are formed by the parameter combi- 
nations for which the discriminant D{aA^ oib] 7a, 7b) is positive. The corresponding opinion 
dynamics then have a unique equilibrium, which (for the parameter combinations in Fig- 
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ure |7|) is approached monotonically. Figure |9| shows a number of graphs fs^aA^aB^A^iB 
parameter combinations with a positive discriminant. The Figure indicates that for smaU 
values of 7^ and large values of 7^ opinion A dominates in equilibrium, and that the domin- 
ion shift towards the alternative opinion if the fraction of contrarians among the A floaters 
increases and that among the B floaters decreases. 

[Figure 9 about here.] 

For given parameters a a and a^, crossing the boundary of the white area in any direction 
away from the lower left corner leads to the occurrence of a saddle-node bifurcation, now 
however outside the domain [a^, 1 — as] of the functions fs^aA^aBHA^iB (maintaining an 
attracting equilibrium in the domain). Therefore in the black region thus entered, the opinion 
dynamics also is characterized by a unique attracting equilibrium. Proceeding towards the 
upper right corner, the line of parameter combinations (7a,7b) satisfying 7^ + 7^ = 1 is 
crossed. On this line the discriminant i^(o^A, 7a, 7b) is singular, and the corresponding 



opinion dynamics have been analyzed in |4.2[ 1. 

In the black area in the upper right corner, for equal and sufliciently small values a a = 
aB = and sufliciently large and equal values 7a = 7b = 7 it has been derived earlier that 



faA^aBHAHB ^ unique unstable equilibrium j5 = 0.5, which causes the convergence of the 
dynamics towards an attracting periodic orbit of period 2. Neither of the two opinions then 
achieves the deflnite majority. The values a and 7 for which this occurs have been derived in 
4.1 , and are represented in Figures [T] and [M] by white line segments in the upper right corners. 



For these parameter combinations the discriminant of the equation fs^aA^aBnA^isip) ~P = 
is negative and thus has three diflFerent solutions, of which two are situated outside the 
domain [a^, 1 — as]- Continuity arguments imply that this behavior will be maintained for 



parameter combinations sufl&ciently close to these line segments. Figure [TO] illustrates this. 
If the parameter combinations are sufliciently far removed from these line segments but 7^ 
and jb are still relatively large (i.e., for given a a and a^, in the upper right corner), the 
dynamics will converge alternately to a unique equilibrium. I.e., by moving away from the 
manifold determined by the constraints aA = c^b and ja = 1b with large values 7^ and 7^, 
a flip bifurcation occurs in which the attracting periodic 2 orbit collapses to an attracting 
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equilibrium point. This is illustrated by Figure [TT] On the attractor the dominion shifts 
towards opinion B with increasing 7^ and decreasing 7^. 



[Figure 10 about here.] 



Figure 10: Four panes of graphs of functions f^-aA,cxB\iAnB relatively small values a a and 
a^, and with 7^ and 7^ satisfying 7a + 7b = 1-8. The values of 7^ and 7^ are indicated 
below each of the four panes. The color coding of the graphs is as follows: red corresponds 
to 7a = 0.8, for green 7^ = 0.85, for blue 7^ = 0.9, for purple 7^ = 0.95, and for magenta 
7a = 1. I.e., for given (aA,<^B), we traverse the line 7a + 7b = 1-8 from its upper left 
point on the 7b = 1 line to its lower right point on the 7a = 1 line. Above each pane the 
densities for opinion A are presented, as obtained by the corresponding opinion dynamics 



f3;aA,aB;7A,7B^ with initial density p = 0.5. 

[Figure 11 about here.] 



Figure 11: Four panes of graphs of functions f'^-aA^o-B-nA^B relatively small values q^a and 
a^, and with 7a and 7^ satisfying 7a + 7b = 1-7. The values of 7a and 7^ are indicated 
below each of the four panes. The color coding of the graphs is as follows: red corresponds 
to 7a = 0.7, for green 7a = 0.8, for blue 7a = 0.85, for purple 7a = 0.9, and for magenta 
7a = 1. I.e., for given (q^a,<^b), we traverse the line 7a + 7b = 1-7 from its upper left 
point on the 7b = 1 line to its lower right point on the 7a = 1 line. Above each pane the 
densities for opinion A are presented, as obtained by the corresponding opinion dynamics 



fs^aA^aB-HAHB^ with initial density p = 0.5. 



We end our discussion by presenting some additional opinion dynamics f^-aA,cxB]iA,iB 
parameter combinations from both the regions with positive and negative discriminant. We 
consider cases with (aA, olb) = (0.1, 0.4) and with (q^a, <^b) = (0.5, 0.3), i.e., a case with small 
a A and intermediate a^, and one with both a a and as intermediate. The corresponding 



graphs of f^-aA,cxB\iAnB represented in Figures 12 and 13, for several values of 7a and 7^. 
All cases allow for a unique attracting equilibrium. High values of both 7a and 7^ lead to 
alternating convergence. Furthermore, a decrease in the fraction of contrarians among the 
floaters of an opinion increases the density of this opinion in equilibrium. 
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[Figure 12 about here.] 



Figure 12: The left column shows four panes of graphs of functions /3;o.i,o.4;7a,7s5 for different 
combinations of 7^ and 7^. The color coding for the different rows is given in the following 
table, with in the second to fifth column the values of 7^: 





red 


green 


blue 


purple 


magenta 


1st (upper) row: 1.8 


0.8 


0.85 


0.9 


0.95 


1 


2nd row: 1.2 


0.2 


0.4 


0.6 


0.8 


1 


3rd row: 0.8 





0.2 


0.4 


0.6 


0.8 


4th row: 0.2 





0.05 


0.1 


0.15 


0.2 



The right column of the Figure shows the densities of opinion A as generated by the corre- 
sponding opinion dynamics for initial value p = 0.5. 

[Figure 13 about here.] 



Figure 13: The left column shows four panes of graphs of functions /3;o.5,o.3;7a,7s5 different 
combinations of 7^ and 7^. The color coding for the different rows is as follows, with in the 
second to fifth columns the values of 7^: 



7a + 7b 


red 


green 


blue 


purple 


magenta 


1st (upper) row: 1.8 


0.8 


0.85 


0.9 


0.95 


1 


2nd row: 1.2 


0.2 


0.4 


0.6 


0.8 


1 


3rd row: 0.8 





0.2 


0.4 


0.6 


0.8 


4th row: 0.2 





0.05 


0.1 


0.15 


0.2 



The right column of the Figure shows the densities of opinion A as generated by the corre- 
sponding opinion dynamics for initial value p = 0.5. 



5 Conclusions 

The results presented re-establish those derived in [2], which concerned communities of non- 
contrarian and contrarian floaters, and [4j, which studied the combined effects of inflexibles 
and non-contrarian floaters. The distinctive patterns of opinion dynamics are not character- 
ized by complete quantitative detail. Rather, the results intend to point to possible outcomes 
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of opinion dynamics. We conclude that various kinds of dynamics may occur. In case the 
local majority rule followed by the contrarian changes are applied for group sizes L = 1 or 
2, the opinion dynamics generically converges to a unique equilibrium. In case the sum of 
fractions of contrarians for the two opinions is larger than 2, the equilibrium generically is 
approached alternately, otherwise the dynamics generically shows a monotone approach. For 
an opinion to obtain the majority in equilibrium, it is required that this opinion is supported 
by a sufficiently large density of inflexibles in combination with a sufficiently small fraction 
of contrarians, as expressed by condition ([5]). 

Group size L = 3 allows for additional outcomes for the opinion dynamics. For sufficiently 
small densities of inffexibles for both opinions, and in addition sufficiently small fractions of 
contrarians among the ffoaters for the two opinions, the dynamics allows for two attracting 
equilibria, that differ in which opinion has the majority. The opinion that eventually will 
achieve the majority thus is determined by the initial condition, and an opinion that has 
a fraction of contrarians that is sufficiently smaller than that of the alternative opinion, 
may achieve the majority in equilibrium, although initially it may be present as a minority. 
For small values of inffexibles in combination with sufficiently large fractions of contrarians 
among the ffoaters, the generated opinion dynamics causes alternating convergence to a pe- 



riod 2 stable orbit (Figure 10). An increase of the densities of inffexibles or a slight lowering 
of at least one of the fractions of contrarians causes the collapse of the attracting periodic or- 



bit into an equilibrium, but maintains the alternating behavior (Figure 11). For a relatively 
large collection of parameter combinations the dynamics ends up on a unique attracting 
equilibrium, which is approached either monotonically or alternating (see e.g. Figures M 



and 13). An increase in the fraction of contrarians among the ffoaters of an opinion leads to 
a decrease of the density of that opinion in equilibrium. Thus, for an opinion to achieve the 
majority in equilibrium, a small fraction of contrarians among its ffoaters is favorable. 
In [2], the "hung elections" outcome in recent votings has been attributed to the interplay of 
non-contrarian and contrarian ffoaters. Likewise, the present paper may shed a light on the 
(dis) appearance of alternating opinion dynamics. The recent wins by Barack Obama and 
Hillary Clinton in two consecutive pre-elections for the Democratic presidential nomination 
2008 in the United States may point to the beginning of an alternating series of wins and 
losses for both candidates. Whether or not this is the case remains to be seen at the moment 
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of writing, but if it does it may point to considerable fractions of contrarians among the 
floaters on both sides, in the states that show alternating election outcomes. In case the 
alternating pattern converges to a stable period 2 cycle, the uncertainty who will win the 
nomination will linger on until the final decisive event. (Note that since the outcome of an 
election in an alternating environment depends on the moment the election actually takes 
place, it may happen that in subsequent state elections the same winner occurs. This is 
however no indication of sustained major support. Furthermore, in a sequence of alternating 
environments, a large number of subsequent wins for the same candidate seems unlikely.) 
If however the alternating changes are converging to an equilibrium, one of the candidates 
eventually will reach a decisive majority. Due to the sensitivity of politics for influences, a 
change in parameter values may easily occur, either with respect to the densities of inflexi- 
bles or to the fractions of contrarians. This may result in a switch from the one alternating 
pattern into the other one, or even into monotone convergence towards an equilibrium. Al- 
though our framework does not map unequivocally to real communities, we think it may 
hint at possible explanations of outcomes of opinion dynamics. 

In forthcoming papers we plan to continue the study of opinion dynamics, by focusing on 
communities in which more than two opinions are being supported. Also mathematical 
aspects of opinion dynamics for group sizes larger than 3 will be studied. 
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6 Appendices 



Subsections |6.H |6.2| and |6.3| present tables for groups of sizes L = 1, 2 and 3 from which 
the density of the A opinion is derived after apphcation of the l.m.r. and the switch by the 
contrarians, given an initial density p for the A opinion. Each table consists of four columns, 
of which the first four are separated by arrows. The first column gives the possible group 
compositions in terms of inflexibles and non-contrarian and/or contrarian fioaters, the second 
column gives the eflFect of the application of the l.m.r. for the group compositions given in 
the first column. An application is indicated by a horizontal arrow, whose first appearance 
in a table is indexed by "l.m.r."; at other places in the tables this index is omitted. The 
third column then gives the effect of the switches by the contrarians if applicable, where it 
is understood that a contrarian switches into a fioater of the alternative opinion. The final 
column gives the contributions of the effect of the l.m.r. and the presence of the contrarians to 
the density of the A opinion, weighed with the probability of the original group composition 
in the ensemble of all possible groups of fixed size, given the densities a a and of the 
inflexibles for both opinions, the fractions 7^ and 7^ of contrarians among the floaters of 
the A and B opinion, respectively, and the densities p — aA for the A floaters and 1 — as —p 
for the B floaters. The total sum of these contributions yields fL-aA^dB-^iA^Bip)' 
In each table the following notation is being used: 

Ai inflexible of the A opinion, 

Af floater of the A opinion, 

Anc non-contrarian floater of the A opinion, 

Ac contrarian floater of the A opinion, 

Af^^ floater of the A opinion coming from a B non-contrarian floater after application 
of the l.m.r., 

Af^ floater of the A opinion coming from a B contrarian floater after application 
of the l.m.r., 

and 
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Bi inflexible of the B opinion, 

Bf floater of the B opinion, 

Bnc non-contrarian floater of the B opinion. 

Be contrarian floater of the B opinion, 

Bf^^ floater of the B opinion coming from a A non-contrarian floater after apphcation 
of the l.m.r., 

Bf^ floater of the B opinion coming from a A contrarian floater after apphcation 
of the l.m.r.. 
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6.1 Table for group size 1 



A A- 



A A 

— > 
— > 

-^nc — > 
Be — > 



aA 

{l--fA){p-aA) 
Bf 



Af iBi^-O^B-p) 
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6.2 Table for group size 2 



A A 




A A 






2a — ^a){p — o^a) 


A . A 




A . A 




A . P n 








AiBi 






olaolb 




— > 


5 Br^e 






-7s)(1 -as-p) 


Be 




Be 




Ai,Af 


2^^75(1 -as-p) 


A A 




A A 






U - lA) [P - OiA) 


A A 




A A 




A-nc^ Bf 


7a(1 -7a)(p-Q^a)^ 


^nci Bi 


> 


^nci Bi 






<^s(l -7a){p-c^a) 


^no Bfie 


— > 


^no Bfie 






(1 -7a)(1 -7s)(;?-Q^a)(1 


A R 




A R 




A-nc^ Af 


2(1 - jAhsip - - o^i 




— > 




— > 


Bf.Bf 





Bi 




Bi 




Bf,B, 





Bfie 




Bfie 




Bf , Bne 





Be 




Be 




Bf,Af 


1a1b{p - ola){1 -as-p) 


Bi, Bi 




Bi, Bi 









Bi , Byie 




Bi, Bjie 









Bi, Be 




Bi, Be 




Bi,Af 


olb1b{'^ - olb -p) 


Bno Bne 




Bno Bne 









Bjio Be 




Bjio Be 




Bnc^ Af 


7b (1 -7b)(1 - olb -pY 


Be, Be 




Be, Be 






as -pf 
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Table for group size 3 







Ai 


A-i, Ai 




a\ 


Ai 


A A 


Ai 


A A 




?>a\{l--fA){p- ola) 


Ai 




_ Ai 


Ai, Ae 


— ^ Ai, Ai, Bf 




A, 


Ai,Bi 


Ai 


Ai, Bi 




2a\aB 


Ai 




Ai 


Ai,Af 




Za\{l -7b)(1 -OLB-p) 


Ai 


Be 


Ai 


Ai^Af^ 


— *■ Ai, Ai, Bf 


2Q!17b(1 - as -p) 


Ai 


A A 


Ai 


A A 




^aA{l-lA?{p-aA? 


Ai 




, Ai 


Anc-) Ac 


— ^ Ai, Anc^ Bf 


4aA7A(l - 1a)(p - OLA) 


Ai 




Ai 


Anc^ Bi 




AaAO^BK^ - lAAP - OLA) 


Ai 




Ai 


Anc^ Af 




6q^a(1 -7a)(1 -1b)\V- ola)^ 

(1 -OLB-V) 


Ai 




_ Ai 


Anc, Af^ 


— *■ Ai,Ane^Bf 


Aola1b(X - 1a){v - OLa)(\ - OLE 


Ai 




Ai 


Ae, Ae 


Ai,Bf,Bf 


OLAl\{V-OLAf 


A, 




Ai 


Ac, Bi 


— ^ Ai,Bf,Bi 


2aAOLB^A{p- ola) 


Ai 


A B 


Ai 


Ae,Af 


Ai,Bf,Af 


4aA7A(l - 7b) b - ola){^ - OLB 


Ai 


A-c, Be 


Ai 


Ac^Af^ 


Ai,Bf,Bf 


20LAlAlBKP - OLA)\X- OLB -V) 


Ai 


Bi, Bi 


Ai 


Bi, Bi 




OLAO?B 


A, 


Bi , BjiQ 


Ai 


Bi, Bnc 




2aAOiB{^ -7b)(1 -olb-v) 


Ai 


Bi, Be 


Ai 


Bi, Be 


— ^ Ai , Bi , Af 


^olaolb1b(X - olb-v) 


Ai 


Bnct Bfic 


Ai 


Bnc J Bnc 




ceA(l -7s)^(l -ces - V? 


Ai 


Bnc-) Be 


_ Ai 


Bnc-) Be 


— ^ Ai,Bnc^Af 


Aola^b(X -7s)(1 - OLB -vf 


Ai 


Be, Be 


_ Ai 


Be, Be 


Ai,Af,Af 


3q^a7|(1 - OLB -V? 
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^nc^ ^nc^ ^nc — > ^nc^ ^nc^ ^nc (1 lfA)^{p (^a)^ 



Anc,Anc,Bi _^ Anc,Anc,Bi 2^^ (1 - 7^)^ (p - Q^a)^ 

Anc,Anc,Bnc A^c^Anc^Af 3(1 - 7a)^ (1 - 7b) " O^a)^ X 

(1 -as-p) 

Anc^Anc^Bc Anc^Anc^Af^ _^ A^c, A^c, B f 27^ (1 - 7^)^ (p - Q^a)^ X 

(1 -as-p) 



Anc^Ac.Ac Anc^Ac.Ac Anc^Bf.Bf {I - 1a)i\{p - ^a)^ 

Bi ^ Anc,Bf,Bi 2q^s7a(1 - 7A)b - Q^a) 

Anc^Ac.Bnc Anc,Ac,Af _^ A^c^Bf^Af 47a(1 -7a)(1 -7s)(p-ceA)^x 

(1 -as-p) 

Anc^Ac.Bc Anc^Ac.Af^ _^ Anc^Bf.Bf 2jAlB{l-lA){p-aA)'^X 

(1 - as -p) 

Anc,Bi,Bi — ^ Bf,Bi,Bi 
A-nc •) Bi , BjiQ — ^ Bf , Bi , BjiQ 

Anc^Bi.Bc Bf,Bi,Bc Bf,Bi,Af 2aBjB{l - 1a){p - (^a)x 

(1 -as-p) 

A-nc 1 Bfic 1 Bfic — > Bf , BjiQ , Byic 

Anc^Bnc^Bc Bf,Bnc,Bc Bf,Bnc,Af 27b (1 - 7a) (1 - 7b) (p " Q^a) X 

{I -as -pf 

Anc. Be, Be Bf, Be, Be Bf, Af, Af 2(1 - jAhUP " «a)(1 - c^b - pf 

Ae,Ae,Ae Ae, Ae, Ae Bf,Bf,Bf 

Ae,Ac,Bi _^ Ac,Ac,Bi _^ Bf,Bf,Bi 

Ae,Ae,Bnc Ae,Ae,Af _^ Bf,Bf,Af {I - -TBhliP - ^a)'^ {1 - aB - p) 

Ae,Ae,Be Ae, Ae, Af^ _^ Bf,Bf,Bf 
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Ac,Bi,Bi _^ Bf^,Bi,Bi _^ Af,Bi,Bi o^\ia{p-ola) 

Ac,Bi,Bnc Bf^,Bi,Bnc Af,Bi,Bnc 2as7A(l -7s)(p-Q^a)x 

(1 -as-p) 

Ac,Bi,Bc Bf^,Bi,Bc Af,Bi,Af AaB'jAjBip - OiA){'^ - o^b - p) 



Ac,Bnc,Bnc — > Bf^,Bnc,Bnc — > Af,Bnc,Bnc 7A (1 " 7^)^ " X 

{1-aB-p)^ 

Ac,Bnc:Bc Bf^,Bnc:Bc Af^B^c^Af 47^7^ (1 " 7s) " Q^a) X 

{1-aB-p)^ 



Ac, Be, Be Bf^,Bc,Bc Af,Af,Af SjajUp - ^a){1 - c^b - p)'^ 

Bi,Bi, Bi — ^ Bi,Bi,Bi 

Bi , Bi , B-YiQ — ^ Bi , Bi , -677,0 

Bi,Bi,Be Bi,Bi,B^ _ Bi,Bi,Af al^feil - as - p) 

Bi-, Brjq^Qj BjiQ — ^ Bi, BjiQj BjiQ 

Bi , BjiQ , Bq — ^ Bi , Byic •) Bq — ^ Bi , B^ic 1 Af 2q;b7b(1 - 7b)(1 - as -p) 

Bi,Be,Be _^ Bi,Be,Be _^ Bi,Af,Af 2aBll{l - ae - pf 

Bno Bnc: -^nc — > B^ic, -^nc? -^nc 

Bnc: Bnc, Be — > ^nc, ^nc, — > ^nc, ^nc, ^/ 7^(1 — 7^)^(1 — OLB — p)^ 

Bnc:Bc,Bc Bnc:Bc,Bc B^c^Af^Af 2(1 - 7^)7! (1 " " 2^)^ 

Be, Be, Be Be, Be, Be Af,Af,Af Isi'^ - (^b - p)^ 
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6.4 L = 3: analysis for the fully symmetric case aA = OiB sind 

The derivative (with respect to p) fs;a,a;j,j the function fs^a^a]^,^ as given by expression ^ 
in the equihbrium p = 0.5 equals 0.5(3 — 4a) (1 — 27). There are two additional equilibria 

^ ^ 1 - 27 ± V(l - 27)(-2 + (3 - 4a)(l - 27)) ^ , ^ . 

P3;a,a;7,7 = ^ 7^ G [a, 1 - Q^J (13) 

1 — Z7 

if and only if < 7 < | and < ^zf^ — 4a. If the two additional equilibria exist they are 
symmetrically positioned on opposite sides of 0.5, and asymptotically stable; the equilibrium 
0.5 then is unstable, with fs-a,a;-f,-f{^-^) ^ 1- 

For a and 7 such that | < 7 < 1 and < — 4a, the equilibrium 0.5 also is unstable, 

> 

with /3;c,a;7,7(0.5) < —1. In this case the dynamics f3-a,a;j,j has two asymptotically stable 
periodic points pl.^,-^ of minimal period 2, symmetrically positioned with respect to 0.5: 

n 1 - 27 ± V(l - 27)(2 + (3 - 4a)(l - 27)) ^ , ^ , 

%;a,a;7,7 = ^'^ ^^7^ G [a, 1 - aj. (14) 

6.5 L = 3: analysis of the general case 

> 

The possible equilibria for h-aA^as-^iA^iB (i^ Wa, 1-Q^b]) follow from solving h-aA.aB-.iA.iB (p) = 
under the restriction that p G [a^, 1 — a^]. We distinguish several cases. 

1. (iaHb) = (0.5,0.5): expression ^ equals /3;aA,«s;0.5,o.5(l^) = 0.5(1 + a^ - o^b), and 
allows for a unique stable equilibrium p = 0.5(1 + a^ — a^), on which opinion A has 
the majority if and only if a^ > a^. 

2. (7a, 7b) 7^ (0.5,0.5), 7a + 7b = 1^ the function fs^aA^aBHA^iB quadratic in p. The 
discriminant 

D{aA,aB]7AnB) for the equation h-aA^cxB-^iAnsip) - P = equals 4a| + 47^(1 - a^ - 
3a|) — 47^(1 — 2(a^ + ol\)). The opinion dynamics /3;a^,a5;7^,i-7^ has a unique equi- 
librium 

1 - 7a + (1 - 27A)a^ - V(l - 7^)(1 - 27A)a| + 7a(1 - 7a) - 7a(1 - 27^)0^ 

(l-27A)(l + aA + aB) 
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in the interval [a^, 1 — as]- The derivative in the equihbrium equals 

1 - 2v/7a(1 - 7a) + (1 - 7a)(1 - 27a)«| - 7^(1 " 27A)«i. 

3. 7a + 7b 7^ 1- The expression fs^aA^aBnA^isip) — P = for determining the equilibria 
now is 

fS]aA:(^B',jA,lB{p) ~ P ~ 

- 7a) + (1 - c^Bhs - (l + 2q^a(1 - 27a) - 7a + 7^)^^+ (15) 
(3 + aA{l - 27a) - asil - 27^) - 47^ - 27^)^? _ 3(^1 _ _ ^^^^3 _ 



Its discriminant is L^(aA,Q^B; 7a, 7b) = (|<?i(<^a, Q^b; 7a, 7b)) + (|<?2(<^a, Q^b; 7a, 7b)) 
with 

Co(QfA, Q^b; 7a, 7b) = Q^a(1 - 7a) + (1 - Q^b)7b, 

Ci(aA,o^B;7A,7B) = -(1 + 2aA(l - 27a) - 7a + 7b), 

C2(q^a, c^b] 7a, 7b) = 3 + aA{l - 27^) - ^^(1 - 27^) - 47^ - 275, 

C3(q^a,Q^b;7a,7b) = -2(1 - (7a + 7b)), 

and 

2 /C2(q^a,Q^b;7a,7b)\^ 1 C2(o^a,Q^b;7a,7b) ci(aA,Q^B;7A,7B) , 

(7i(aA,Q^B;7A,7B) = ^ ^ \ + 

27 Vc3(aA, c^B] 7a, 7Bj ^ 3 03(0^^, c^b] 7a, 7Bj ^3(0^^, Q^b; 7a, 7Bj 

Co(aA,Q^B;7A,7B) 

C3(aA,Q^B;7A,7B)' 

1 /C2(aA,Q^B;7A,7B)\^ , Ci(aA,Q^B;7A,7B) 
(?2(o^a,Q^b;7a,7b) = -- ^ ^ + 



3 \Cs{aA, o^b; 7a, 7b) ^ ^3(0^^, Q^b; 7a, 7b) 



Figure 14 shows, in three parts, a signplot of the discriminant for this case. In addition 



the results of the analysis for combinations (a, a; 7, 7) is included. 

[Figure 14 about here.] 
[Figure 15 about here.] 
[Figure 16 about here.] 
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Figure 1: Figure a shows the graphs of /i;o.2,o.2;o,o (which coincides with the diagonal) and 
/i;0.2,o.2;0.075,o.05 fuuctious of p ou the interval [0.2,0.8]. By changing the frequencies of 
contrarians from (7^,7^) = (0,0) into (7^,7^) = (0.075,0.05), the collection of neutrally 
stable equilibria (the diagonal) bifurcates into a unique stable equilibrium p — 0.44 on which 
the B opinion has the majority. Figure b shows the diagonal together with the graph of 
/i;0,0;0.2,o.2 0^ [O, l], and the graph of /i;o.i5,o.i75;0.2,o.2 on [0.15,0.825], both as functions of p. 
The dashed lines indicate the boundaries of the interval [0.15, 0.825]. The graphs of the two 
functions almost coincide on this interval. The dynamics generated by these two functions 
have p = 0.5 and p = 0.4875 as their respective stable equilibria. 
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Figure 2: 
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/3;0.1,0.1;7,7 




Figure 3: Graphs of /3;o.i,o.i;7,7 as function of p G [0.1,0.9], with values 7 as indicated at 
each specific graph. In addition, the diagonal and the line 1 — p are drawn. 
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gure 4: The critical curves C3 and C3 in the (a, 7)-parameter space 
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0.5 0.5 

7A = 7A = 0.1 



0.5 

7A = 0.2 



0.5 

lA = 0.3 



0.5 

7A = 0.4 




Figure 5: The different panes, distinguished by different values of 7^, have a a on the 

horizontal axis and as on the vertical one. Each pane shows the line aA + c^b = 1 and 

1 - 27a 1 + 27a 

in addition the dashed line aB — \ c^a of parameter values for which the 

3 - 27a 3 - 27a 

equilibrium value Ps;aA,aBnA^-iA equals 0.5. Below the dashed lines the equilibrium value 
lies above 0.5, i.e., opinion A then obtains the majority. The region q^a + c^b > 1 is not 
involved in the analysis. 
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Figure 6: The different panes show in white the region of parameters (q^a^Q^b) for which 
the equihbrium P2>-aA,oiB\iA^-iA approached monotonically; the black regions indicate pa- 
rameter combinations for which the equihbrium is approached alternately. For 7^ = 0.5, the 
derivative of fs^aA^aB^A^-iA equilibrium equals for all parameter values (QfA,Qf^), 

and the equilibrium is reached in one iteration. The line aA + o^b = 1 refers again to the 
degenerate cases for which the opinion dynamics of density p is restricted to a single density 
p — Parameter combinations (q^a^Q^b) for which aA + c^b > ^ are not involved in the 
analysis. 
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(c^A, c^b) = (0, 0) (q^a, c^b) = (0.1, 0) 



Figure 8: 
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Figure 9: Four panes of graphs of functions fs^aA^aB-^AnB relatively small values a a and 
as^ and with 7^ and 7^ satisfying 7a+7b = 0.7. The values of 7^ and 7^ are indicated below 
each of the four panes. The color coding of the graphs is as follows: red corresponds to 7^ = 0, 
for green 7^ = 0.2, for blue 7^ = 0.35, for purple ja = 0.5, and for magenta ja = 0.7. I.e., 
for given (a^, a^), we traverse the line 7a + 7b = 0.7 from its upper left point on the 7a = 
line to its lower right point on the 7^ = line. The discriminant i^(o^A, 7a, 7b) for the 
exposed parameter values is positive, indicating a unique equilibrium for the corresponding 
opinion dynamics. 
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Figure 10: 
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Figure 11: 
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Figure 12: 
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Figure 13: 
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Figure 14 a 



aA = 0.2 
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as = 0.5 
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Figure 14 b 



aA = 0.5 
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Figure 14: Figures b, and c shows a collection of panes with signplots of the discriminant 

a^; 7^, 7^), for values aA and as as indicated, and 7^ and 7^ for each pane ranging 

between and 1, with 7a + 7b 7^ 1; 7a is on the horizontal axis and 7^ on the vertical 

axis. The white areas represent the parameter combinations with a positive discriminant, 

and in dark regions the discriminant is negative. On the line 7^ + 7b = 1 in each pane the 

discriminant becomes singular. In panes for which aA = c^b holds, on the line 7^ = 7b in 

the black regions (i.e., a negative discriminant) in white the points are indicated for which 

> 

the equilibrium p = 0.5 for fs;aA,aBnA,iB i^ unstable; other points on the lines ja = 1b 
(for a A = oib) indicate parameter combinations for which p = 0.5 is stable (as follows from 
Figure ^ . 
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